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ABSTRACT 

An open loop optimal control algorithm is developed for 
general flexible structures, based on Laplace transform 
methods. A distributed parameter model of the structure is first 
presented, followed by a derivation of the optimal control 
algorithm. The control inputs are expressed in terms of their 
Fourier series expansions, so that a numerical solution can be 
easily obtained. 'Hie algorithm deals directly with the 
transcendental transfer functions from control inputs to outputs 
of interest, and structural deformation penalties, as well as 
penalties on control effort, arc included in the formulation. The 
algorithm is applied to several structures of increasing 
complexity to demonstrate its generality. 


L INTRODUCTION 

The control of large flexible structures has become an 
important issue in recent years, primarily in the aerospace 
industry. 1 As larger structures continue to be deployed in 
space, the effects of control-structure interaction arc becoming 
increasingly important. For example, stringent pointing 
requirements for space-based antennae make it necessary to 
isolate and suppress unwanted structural vibration caused by 
both slewing maneuvers and exogenous disturbances. 
Consequently, it becomes necessary to model structural 
flexibility when developing control laws for these types of 
structures. 

Because disturbances and control forces generally act at 
discrete points on the structure, structural responses tend to 
exhibit wave propagation characteristics. Traditional finite 
element codes are unable to capture the high frequency behavior 
of such structures, due primarily to the spatial discretization 
associated with lumped parameter models. This limitation 
makes it particularly difficult to study the propagation of 
flexural waves within structures, since an extremely fine 
discretization is required to preserve the local wave- like 
characteristics of the disturbances. To overcome this problem, 
this paper develops a distributed parameter, system-based 
model, which deals directly with the governing partial 
differentia] equations that describe the structure. 

Given the continuum model of a flexible structure, there 
remains the issue of identifying control methodologies that take 
advantage of the additional high frequency information 
available therein. Tzafcstas 2 develops a distributed parameter 
analogue of the linear quadratic regulator theory. A distributed 
parameter Riccati equation, expressed in terms of spatial 
differential operators, is presented. Miller, Hall, and von 
Flotow^ develop optimal control laws for power flow at 
structural junctions based on a travelling wave approach. The 
effect of the localized controller is to modify the wave scattering 
matrix at the junction in a way that minimizes the power 
flowing from the junction. MacMartin and 1 tall 4 consider 

*Mcniber, AIAA 


optimal control of power flow in uncertain structures based on 
an Hoo cosl criterion. Closed-loop stability is guaranteed by 
minimizing the maximum power imparted to the structure over 
all frequencies. The optimal distributed control of a rigid 
spacecraft with flexible appendages is discussed by 
Break we 11. ^ 

Skaar 6 presents closed-form open loop optimal control 
solutions for a simple structure. The cost function considered 
has the form: 

h 

J - J { k]ii(l) 2 + k2u(t) 2 } dt (1) 

where u represents the control input, and k\ and k 2 arc 
constants. Terminal and integrated penalties on the structural 
deformations are not permitted. Rather, the terminal constraints 
arc adjoined to the cost function with Lagrange multipliers. 

The exclusion of dcformational penalties makes it possible to 
derive analytical solutions for certain types of maneuvers. 
Otherwise, the optimal control solution can, in general, only be 
obtained by numerical methods. 

Analytical results arc available for only the simplest of 
distributed parameter models, containing very few flexible 
elements. More often, a complex structure, such as a truss 
beam, is replaced by a single equivalent member in the 
continuum model. Such an approximation is usually accurate at 
low frequencies only. For general structures, the structural 
responses must be calculated numerically. The convolution 
integral representation technique developed by Skaar is 
generalized in this paper to handle arbitrary structural 
configurations. 

A review of the continuum modeling approach is 
presented in section 2. The optimal control formulation is 
developed in section 3. Several examples of this method, 
applied to structures of increasing complexity, are then 
presented in section 4. Conclusions and recommendations can 
be found in section 5. 


2 ■ STR UCTURAL MODELING 
2A Modeling of Flexibl e Elements 

I raditional approaches for modeling complex structures 
have relied on finite element modeling techniques. This 
approach idealizes a structure as an assembly of many small 
pieces which arc constrained to move together in a manner 
consistent with the internal clastic behavior of the underlying 
continuum model. These techniques arc powerful and widely 
used. I lowcvcr, they suffer from various modeling 
idealizations which limit the accuracy of behavior predictions, 
particularly for high frequency. 

In order to better model the high frequency behavior of 
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clastic fra me- 1 ike structures, a continuum approach is presented 
in this paper which overcomes the conventional limitations of 
traditional finite clement modeling techniques. The continuum 
method (also known as the exact finite element, distributed 
parameter, or dynamic stiffness method), deals directly with the 
governing partial differential equations for the individual clastic 
elements to eliminate the explicit time dependence in the 
equations of motion. The Laplace transform is employed to 
convert the governing partial differential equations into ordinary 
differentia] equations in the spatial dimension. For common 
element models (e.g., rods in torsion, Bcmoulli-Euler beams in 
bending), simple analytical solutions to these equations exist. 
The resulting solutions are explicit functions of the generally 
complex frequency parameter, s, which has been introduced 
through the application of the Laplace transform technique. 

As an example, consider the case of a rod in torsion, 
shown in Fig. 1. The governing partial differential equation is 


GJe"(x,t) + mr 2 0(x,t) = x(x,t) 


( 2 ) 


where 0 is the cross sectional angle of twist, GJ is the torsional 
rigidity, m is the mass per unit length, r is the cross sectional 
radius of gyration, ( )’=r)( )/dx, ( )=d( )/0t, and x is the 
distributed torque. Applying the Laplace transform, we obtain 

0 "(x,s) + tn £rj L ~0(x,s) = 0 (3) 

where initial conditions and distributed forcing along the 
clement have been temporarily neglected. The general solution 
follows as: 

0(x,s) = A(s) cos px + B(s) sin px, P 2 = ‘~qJ~ s2 


where the functions A(s) and B(s) depend on the constraints 
imposed at the boundaries of the element. To simplify the 
mathematical developments, a structural state vector is defined 
for the elastic element as follows: 


y(x,s) 



(5) 


where T represents the net torque resultant along the rod. With 
knowledge of the state at one boundary of the element, it is then 
possible to determine the state at any internal location. A 
frequency-dependant spatial transition matrix is used to 
propagate the system state to arbitrary clement locations. Its 
transcendental elements can be thought of as spatially varying 
transfer functions. For this example, the transition matrix can 
be shown to be 


d>(x,s) = 


cos px 
GJp sin px 


gLsinpx 
COS px_ 


(6) 


and 


y(x,s) = d>(x,s) y(0,s) 

Alternatively, when the displacements at the boundaries are 
known, they can be related to the forces. This is accomplished 
by a dynamic stiffness matrix. Its name derives from the 
stiffness matrix associated with the traditional Finite clement 
method. For the rod in torsion, this matrix is given by 


K(s) 


GJg COS pi -1 
sin pi ^ _1 cos pi 


(8) 


where 1 is the length of the rod. 

In addition to torsional rods, the formulation also 
handles Euler beams in bending in two directions and axial 
rods. Timoshenko beams can also be included as continuum 
elements. 

Because the continuum method treats each elastic 
member as a single clement, no spatial discretization is 
required. This is in marked contrast with traditional finite 
clement methods, where each element must be lumped into 
several segments. As a result, the continuum model is valid at 
all frequencies, insofar as the partial differential equation 
represents the actual physical structure. In contrast, the spatial 
discretization associated with the finite element approach creates 
a computational burden for even the simplest of structures (each 
flexible element is typically broken down into more than ten 
segments). This significant reduction in the number of discrete 
modeling elements required makes the continuum method more 
attractive from a computational point of view. 

Another advantage to the frequency domain modeling 
approach is the ease with which damping is incorporated into 
the structural model. For internal damping (where energy is 
dissipated as heat within the structural elements) the static 
bending, axial or torsional stiffness is replaced by a complex 
valued function of the complex frequency. The functional 
relationship depends on the type of damping modeled. For 
example, a fractional derivative damping model scales the static 
stiffness by the square root of the complex frequency. 7 This 
type of damping model is extremely difficult to implement in 
time domain formulations, and requires a knowledge of the 
entire past history of the deformation of the structural element. 
For external damping (where structural energy is dissipated to 
the surroundings) the mass per unit length is replaced by a 
frequency-dependant parameter. 7 


2.2 Assemb ly of Elements 

The assembly of flexible elements into a complex 
frame-like structure is accomplished using the method of 
local/global coordinates, 7 which is implemented in most finite 
element software. The structure to be modeled is divided into a 
set of flexible elements and a set of rigid joints, which attach to 
any number of flexible elements at their respective boundaries, 
as shown in Fig. 2. External forces arc applied at the joints 
only, but the deformation of the structure is available at all 
points. (The case of a concentrated force applied within a 
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Fig. 1 : Example of a distributed parameter element: a rod in torsion. 



Fig. 2: Generic frame- like structure. 


flexible element is treated by breaking the member into two 
continuous elements connected by a massless joint. Modeling 
distributed forcing is somewhat more difficult.) Massive joints 
are modeled by dynamic stiffness matrices as well, so that their 
contribution to the structural response at all frequencies is 
retained. 

The topology of the structure is given by a connectivity 
matrix, which relates the local displacements of the structural 
elements (both flexible and rigid) to a set of global 
displacements which uniquely describe the location and 
orientation of all joints in the structure. The applied forces at 
the joints arc defined in a dual manner, so that for every global 
displacement there exists a global forcing input at the same 
point and in the same direction. It is then a simple matter to 
compute tiie dynamic stiffness matrix for the structure as a 
whole at any given complex frequency. The individual 
stiffness matrices are first arranged in a large block diagonal 
matrix. This matrix is then post- and pre-multiplicd by the 
connectivity matrix and its transpose, respectively, resulting in 
the system dynamic stiffness matrix. 

Because the dynamic stiffness matrix is transcendental 
in nature, computing natural frequencies is not a simple matter 
of solving an eigenvalue problem, as is the ease for the finite 
element approach. Rather, the stiffness matrix must be 
computed at many frequencies in order to gradually converge 
on each modal frequency. However, a powerful algorithm is 
available which rapidly converges on these eigenfrcqucncics. 8 
The algorithm works for undamped systems only, and 
additional root- searching algorithms must be employed when 
structural damping is modeled. 7 

The system dynamic stiffness matrix can be inverted, if 
desired, yielding the dynamic flexibility matrix for the 
structure. Its elements can be thought of as transfer functions 
from the joint forces to the joint displacements or as admittance 
functions. This flexibility matrix, in conjunction with the 
flexible element spatial transition matrices, enables the 
straightforward calculation of the transfer function from any 
joint force to any point on the structure. 


2.3 Inverse Laplace Transform 

In the frequency domain, the responses at various 
locations within a linear elastic structure to multiple control 
inputs are determined by exploiting the principle of 
superposition. The set of responses, yi(s), are expressed as 

m 

y»(s> = I gij(s) Uj(s), i = 1 r (9) 

j=l 

where gij(s) is the transfer function from the j’th input, uj(s), to 
the i’th output, yj(s). These equations can be expressed in 
more compact form using matrix notation: 

y(s) = G(s)u(s) (10) 

The matrix G(s) is the dynamic flexibility matrix for the given 
structure, or some partition of it, depending on the inputs and 
outputs considered. 

The time response vector corresponding to y(s) is 
available via the inverse Laplace transform, given by 

a+j~ 

y(0 = 2 b J y (s) e st ds (11) 

J a-j«» 

where the integration path is known as the Bromwich contour. 9 
For finite dimensional systems, a residual expansion is used in 
lieu of Eq. (1 1) to compute the time responses analytically. 


The same can be done for distributed parameter systems, except 
that the expansion has an infinite number of terms and must 
therefore be truncated at some point. However, greater 
numerical accuracy is possible by working with Eq. (11) 
directly. The frequency domain response is tabulated for 
values of complex frequency equally spaced along the 
Bromwich contour, and a numerical procedure converts this 
data into a response history evaluated at equal spaces in time. 

Because this approach utilizes frequency domain 
representations of the control inputs, it circumvents the 
computationally expensive calculation of convolution integrals. 
Furthermore, signals that cannot be represented by finite 
dimensional state space models are easily handled in the 
frequency domain. For example, the implementation of a time 
delay simply requires multiplication of the frequency domain 
data by a suitable exponential of the complex frequency before 
the inverse Laplace transform algorithm is invoked. 
Implementing such a time delay on a modal basis requires a 
truncated series expansion of the complex exponential (such as 
a Padfe approximation), with many terms needed to obtain an 
accurate representation. 


3 OPTIMAL CONTROL FORMUL ATION 

With the continuum modeling approach described 
above, it is possible to recast a class of optimal control 
problems into a convenient form, from which optimal control 
trajectories are easily calculated. This form is applicable to a 
completely general frame-like structure (although applying this 
method to structures containing plates and membranes is the 
subject of current research 1 !), with multiple control inputs and 
multiple outputs. The class of problems discussed here are 
fixed-time, linear quadratic, open loop control problems with 
penalties on control effort, position and velocity of various 
output points on the structure, and structural deformation. 

Thus, the cost functional has the form 

J = ly(tr)- y<j] T ii[y(tf)- ydl 

if 

+ J ( ly(tf) - ydl T Qly0r)- yaJ + u(t) T Ru(t)}dt (12) 
0 

where II, Q and R are weighting matrices, and y t | is the vector 
of desired output values. It should be noted that y t i represents 
the physical output variables of interest, and is not related to the 
outputs of some state variable representation of the system. 

Traditionally, the dynamics of the system are adjoined 
to this functional via a costate vector as differential equation 
constraints. However, because the structural transfer functions 
are transcendental and infinite dimensional, a finite dimensional 
costate vector cannot be defined. Tzafestas 2 succeeded in 
identifying a distributed parameter optimal control solution 
which incorporates an infinite dimensional costate. This 
solution represents the distributed parameter analogue of the 
Riccati differential equation for finite dimensional systems. 
However, the method is not immediately applicable to complex 
structures, where more than one partial differential equation is 
involved. Even the case of a single beam in bending presents 
considerable difficulty. 12 

One alternative to adjoining an infinite dimensional 
costatc is modal truncation. The high frequency modes of the 
structure are simply ignored, and the dynamics of the structure 
is approximated with a finite dimensional state space 
realization. However, in order to take advantage of the 
"exactness" of the continuum modeling approach, it seems 
appropriate to avoid modal truncation altogether. Instead, we 
express each response as the convolution of impulse responses 
with control inputs. This yields for the cost function 
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“if it r i r 

J= j G(t r -x)u(x)dt- y ( i II J G(tf-x)u(t)dx - y d 

't 

+ [ [ [ J G(t-x)u(x)dx- y d Q | (J(t-x)u(x)dT -y d 


+ u(i) Ru(t)j dt (13) 

where G(t) is the matrix of impulse responses from each 
control input to each output, as defined in section 2.3. These 
convolution integrals need not be computed directly, as they arc 
the inverse transforms of the appropriate transfer functions 
multiplied by the associated control inputs. Taking variations 
in u, we arc left with an integral equation which, except fora 
very small number of special cases, is difficult or impossible to 
solve in closed form. It is therefore necessary to express the 
control inputs as weighted sums of an appropriate set of basis 
functions which span the function space of allowable control 
inputs. For fixed time problems, the Fourier scries is a good 
candidate, leading to 


f(0 T c m 


= m e, f(o = 


sin rct/tf 

COS Kt/tf 

sin nrciAf 
-cos rntt/tf 


The cost functional now depends only upon a constant vector, 
c, representing the coefficients in the Fourier series expansions 
of the control inputs. The resulting cost functional is quadratic 
in the coefficient vector, and the minimization problem is 
straightforward, yielding: 


c = A' Byu 


A = Y(t f ) T IIY(tf) + J ( Y(t) T Q Y(t) + F(t) T RF(t) } dt (17) 
0 


B = Y(tf) T II + J { Y (t) T Q } dt 
0 


Y(t) = J G(t-x)F(t)dx 0 ' 

Once again, the convolution of impulse responses with basis 
function inputs can be calculated via the inverse Laplace 
transform: 


Y(t) = L- | (C.(s)F(s)) 


Furthermore, these basis responses can be computed a priori, 
provided that the forcing locations and structural deformation 
penalty locations are known in advance. This makes it possible 
to try a large number of cost functionals without repeatedly 
calculating the responses to basis inputs. It should be 
mentioned that a large amount of memory is required to store 
this data. 

A unique advantage of this approach is that it readily 
accommodates penalties in higher derivatives of both control 
effort and physical deformation. In the frequency domain, 
differentiation merely requires multiplication of the data by the 
Laplace transform variable. The inverse transformation will 
then produce the derivative of the original signal. Higher order 
derivatives are obtained by multiplying by higher powers of the 
complex frequency. Incorporating higher derivative penalties 
in the traditional optimal control formulation is considerably 
more difficult. 

Special constraints on the control histories arc treated by 
adjoining the constraints via Lagrange multipliers. An example 
is the requirement that the controls be continuous at the 
beginning and end of the maneuver. This implies that the 
controls go to zero at the initial and final times. Two Lagrange 
multipliers arc therefore introduced for each control input. 
Performing the minimization, we obtain 

-c -I A F(0) F(tf)r l rBy d - 
>-1 = F(0) T 0 0 0 (21) 

_*2J LF(t f ) T 0 0 J [o 

where and X 2 arc Lagrange multiplier vectors corresponding 
to the control constraints at the initial and final times, 
respectively. 

It should be noted that the only approximation in the 
entire development involves expressing the control inputs in 
terms of the basis functions. The dynamics of the entire 
structure is accounted for, since the impulse responses are exact 
(insofar as the original equations represent physical reality). 
Also, the structural deformations are assumed to be small, so 
that linearization does not introduce significant errors. As a 
result, large angle slew maneuvers are not included in this class 
of problems. It is possible, however, to express structural 
deformations with respect to a nominal condition during a large 
angle slew, and then linearize about that reference. 

It is important to note the difference between the 
approximations made in the continuum approach and those 
made in control systems based on finite clement models. 
Typically, a finite element model is used to determine a 
truncated state space realization of the system. The control 
system is then designed using standard methodologies, such as 
linear quadratic regulator (LQR) or theory. The modeling 
error associated with modal truncation must then be considered 
in assessing the performance of the system. In contrast to this, 
the continuum approach retains all modes of the system and 
therefore has no modeling error associated with it, provided 
that the original partial differential equations represent the 
physical system exactly. The only limitation to this approach is 
therefore the number of basis functions chosen to represent the 
control signals. Using the Fourier series expansion essentially 
places a limit on the bandwidth of the control inputs. If the 
bandwidth of the physical controls are known in advance, the 
number of basis functions used in the calculations can be 
adjusted accordingly. The resulting control history would then 
indeed be optimal for that set of physical controls. 
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4. EXAMPLES 


10.0 


4.1 UiiMd Mass with E ksiMfi Annendaue 

In an earlier analytical study by Skaar 6 , the open loop 
control of a rigid mass with a flexible appendage, shown in 
Pig. 3, was studied. In his work, dcformational penalties were 
not incorporated into the cost function; rather, the terminal 
conditions were adjoined to the cost functional as constraints. 
Skaar derived analytical expressions for impulse responses of 
the simple mass/appendage structure and thus obtained closed 
form optimal control solutions for the structure. His approach 
does not readily generalize for more complex structures. In 
contrast, the formulation presented here readily generalizes for 
realistic complex structures. Nevertheless, the mass/appendage 
structure is used as a first example to validate the optimal 
control formulation. 

The maneuver involves translating the mass a distance 
of 10 meters along the axis of the flexible appendage, bringing 
it to rest with minimal residual energy and post-maneuver drift 
after 20 seconds. The first case places terminal penalties on the 
final position and velocity of the rigid mass and on a point 4/5 
of the length along the flexible appendage. A small penalty is 
also placed on control rate, and 17 basis functions are used to 
approximate the control input. The results, shown in Fig. 4, 
indicate that the terminal conditions are matched, and residual 
energy is minimized. The second case places additional 
terminal penalties on points along the flexible appendage, as 
well as en -route (integrated) penalties on structural deformation 
at these points. As shown in Fig. 5, the magnitude of the 
structural deformation during the maneuver is reduced slightly 
(about 15%), and large deformation occurs over smaller 
intervals of time, resulting in a lower RMS displacement of the 
tip of the appendage over the maneuver interval. I he terminal 
conditions are again matched. In the final case, the member 
stiffness is reduced by a factor of four, so that the primary 
modal frequency of the structure corresponds approximately 


X(t) 


n o 

? 


C 



(b) 



(c) 

Fig. 4: Results of optimal maneuver of mass/flcxiblc appendage system 
(case #1): (a) control force applied to rigid mass, (b) position of rigid mass, 
(c) deformation of tip of flexible appendage with respect to position of 
rigid mass. 


u(t) 



M = 1 kg 


3 


l = 0.5 m 
m = 2 kg/m 
EA = 0.05 Nt 


Fig. 3: Mass/flcxiblc appendage system. 


with the frequency of the first basis function of the control 
input. The results of this case, presented in Fig. 6, indicate that 
the control input has been adjusted so that excitation of the 
primary mode of the structure is suppressed, at the expense of 
increased control effort. Again, the terminal conditions arc 
matched, and residual internal energy is minimal. 



(a) 


4.2 Starfi&]i_CunfiguratiQn 

A typical model for a spacecraft with flexible 
appendages is the starfish configuration, shown in Fig. 7. It 
consists of a rigid hub, to which four flexible arms arc 
attached. At the end of each ami is a rigid mass. Control 
inputs are available at the hub and at the tips of two of the 
flexible arms. 

For this structural model, two maneuvers arc presented. 
The first is an in-plane translation of 0.1 meter in the direction 
of one of the flexible appendages, with a final time of 5 
seconds and all three control inputs available. Terminal 
penalties are placed on the position and velocity of the hub, as 
well as at three locations along the two arms whose axes lie in a 
direction perpendicular to the motion. It is anticipated that the 
control forces will excite bending vibration in these 
perpendicular members due to their inertia. Identical magnitude 







and velocity penalties arc also placed on the three control 
inputs. The results of this maneuver are shown in Fig. 8. The 
forces provided by the thrusters at the two appendages arc 
identical (except for their sign) and differ slightly from the 
control history of the thruster at the hub. This discrepancy 
compensates for the inertia of the flexible arms, bringing the 
system to the final desired state with minimal residual energy. 

The second maneuver involves a small (0.1 radian) 
counter-clockwise in-plane rotation about the hub in addition to 
the translation of the previous maneuver, in order to 
demonstrate the ability of the optimal control approach to 



Time (sec) 
(c) 


Fig. 5: Results of optimal maneuver of mass/llcxiblc appendage system 
(ease #2): (a) control force applied to rigid mass, (b) position of rigid mass, 
(c) deformation of Lip of flexible apjKndagc with respect to position of 
rigid mass. 





Fig. 6: Results of optimal maneuver of mass/flcxiblc appendage system 
with reduced axial stiffness (ease #3): (a) control force applied to rigid 
mass, (b) position of rigid mass, (c) deformation of lip of flexible 
appendage with respect to position of rigid mass. 


handle multiple final conditions. Here, only the lateral forces at 
the tips of two of the appendages arc available. Also, terminal 
position and velocity penalties arc placed on two points along 
all four arms, as well as on the central hub, and the final time 
remains unchanged. The results are presented in Fig. 9. 

Again, the incorporation of state penalties involving structural 
deformation succeeds in minimizing residual internal energy at 
the terminal time. As previously mentioned, the linear control 
solution applied to this type of maneuver is only appropriate for 
small rotations and angular rates. For large rotations, nonlinear 
kinematics must be considered, whereas for large angular rates, 
gyroscopic forces become significant. 
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Fig. 7: Diagram of starfish con figuration. 





thruster #3 (double dotted line), (b) position or hub. (c) relative transverse 
deflection of flexible arm associated with control #2. 





(c) 

Fig, 9: Optimal maneuver of starfish involving translation and rotation: 
(a) Control forces applied to thruster #2 (solid line) and thruster #3 (dolled 
line), (b) translation of central hub, (c) rotation of central hub. 


4..1 STOLE Structure 

The final example is a complex three dimensional 
structure proposed by NASA as a design challenge.^ The 
Spacecraft COntrol Laboratory Experiment, shown in Fig. 10, 
consists of a rigid shuttle and hexagonal truss antenna 
connected by a flexible mast. Previous authors have treated the 
antenna as being rigid. In this paper, however, the flexibility 
of the antenna is considered. Figure 1 1 shows the transfer 
functions from a torque applied to the shuttle along about the 
axis of the mast to various points along the mast and antenna 
for both the rigid and flexible antenna models. A validation of 
this structural model with a traditional finite element code is 
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currently underway. As indicated in the figure, the flexibility 
of the antenna has a considerable effect on the transfer 
functions at higher frequencies. It is therefore appropriate to 
include a dynamic model for the antenna when applying the 
optimal control algorithm. 




(b) 

1° : The SCOLE experiment: (a) Shuttle and flexible mast 

properties, (b) flexible antenna model used in this study. 






Fig- 1 1 : Some transfer functions from control torque applied to 

shuttle along mast axis to inertial accelerations at various points on 
SCOLE structure for a rigid antenna (solid lines) and for a flexible antenna 
(dolled lines): (a) torsion of midpoint of mast, (b) torsion of mast at 
antenna junction, (c) transverse deflection of mast in pitching direction at 
antenna junction, (d) deflection, in plane of antenna, of antenna hub. 


The SCOLE maneuver presented here consists of a 0. 1 
radian rotation about the z-axis of the shuttle. The three 
dimensional model used in the optimal control formulation 
incorporates 52 partial differential equations, which account for 
axial, torsional, and bending vibration in the mast and each of 
the twelve antenna elements. The shuttle is modeled as a 
massive rigid body with six degrees of freedom. Torque 
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controls directed along the z-axis arc placed at either end of the 
mast. Due to the asymmetry of the structure, bending/torsion 
coupling is expected. Consequently, roll and pitch torque 
controls are also located on the shuttle. In the cost functional, 
equal magnitude and rate penalties arc placed on all four 
controls. Large terminal penalties arc applied to the roll, pitch 
and yaw angles of the shuttle, as well as the torsional 
deformation of the mast at its midpoint and at the mast/antenna 
junction. The maneuver lime is ten seconds. 

The results of the SCOLE slew are shown in Fig. 12. 
Because of the conventions used in defining the axes at the 
nodal points, the torques applied to the ends of the mast arc 
opposite in sign. The physical torques are, however, applied in 
the same direction. From the figure, it is clear that, although 
the shuttle has rotated the prescribed amount, there is a small 
amount of residual torsional energy in the stmeture. This 
energy is due primarily to the deformation of the antenna. It is 
expected that additional penalties on structural deformation 
placed at various points on the antenna will significantly reduce 
tli is residual energy. 


5. CONCLUSIONS 

The open loop optimal control algorithm has been 
demonstrated for several structural models. The generality of 
the approach has been exploited in applying the algorithm to a 
complex structure. Furthermore, the ability of this approach to 
handle constraints and derivative penalties has been 
demonstrated. One major issue concerning this open loop 
method is modeling error. Because the approach presented is 
open loop, this formulation makes no guarantees on the 
performance of the actual structure, for which the mathematical 
abstraction is only an approximation. The necessary sensitivity 
analysis and closed loop formulations are topics of current 
research. It may also be possible to incorporate an adaptive 
identification algorithm to adjust the model in such a way as to 
minimize modeling error. 

If the optimal control problem has an infinite time 
horizon, it is possible to use a different set of basis functions to 
converge on a solution. Candidate bases include Legendre and 
Laguerre functions, which arc linear combinations of 
exponentials and other decaying functions J 4 

A method for incorporating distributed penalties into the 
optimal control algorithm is one topic of current research. 
Rather than penalizing certain points along a particular flexible 
element, an integral over the entire element, subject to some 






(d) 

Pig. 1 2: Results of optimal SCOLE maneuver: (a) yaw torque 

applied to shuttle (solid line) and masl/anlcnna junction (dotted line), (b) 
roll torque (solid line) and pitch torque (dolled line) applied to shuttle, (c) 
yaw rotation of shuttle, (d) relative torsional deflection of mast at midpoint 
(solid line) anti at mast/an lenna junction (dolled line). 


weighting function, becomes part of the cost functional. Such 
a capability would make it possible to develop controls that 
minimize the total (kinetic plus potential) energy within the 
stmeture in a continuous sense, rather than penalizing a large 
number of points within each flexible element. 
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Work is currently underway to extend the basic 
modeling approach to continuous plate and membrane 
elements. This capability will permit modeling and control 
formulations for the NASA space station class of large space 
structures. 
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